Estimating purity in terms of correlation functions 
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We prove a rigorous inequality estimating the purity of a reduced density matrix of a composite 
quantum system in terms of cross-correlation of the same state and an arbitrary product state. 
Various immediate applications of our result are proposed, in particular concerning Gaussian wave- 
packet propagation under classically regular dynamics. 
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The autocorrelation function of time evolution, also 
known as survival probability or as fidelity in the con- 
text of echo dynamics [ij, is an important tool in the 
discussion of quantum systems. Indeed the autocorrela- 
tion function contains a vast amount of information, but 
analysis in terms of the more general concept of cross 
correlation functions can bring additional insight. This 
has been recently demonstrated for spectral statistics . 
In this note we shall seek a better understanding of the 
evolution of entanglement in terms of cross correlation 
functions. This is particularly attractive, because entan- 
glement is at the very root of quantum mechanics, and 
we shall therefore test the usefulness of cross correlation 
functions in this context. Note that the first experimental 
test on the value of cross correlations for the analysis of 
spectral statistics was carried out in a microwave exper- 
iment with classical fields rather than in the context 
of quantum mechanics. 

Entanglement indicates to what extent the state under 
consideration can be written as a product of states in two 
subsystems selected due to their interest in the physical 
context. For pure states, entanglement is reflected in the 
properties of the reduced density matrix for any of the 
two subsystems. While it is tempting to use some form 
of entropy to describe this property, this makes analytic 
work very difficult. We shall therefore use the purity of 
a subsystem as a measure of entanglement. Purity is 
defined as the trace of the square of the reduced density 
matrix. The fact, that we use an analytic function, allows 
to obtain the basic inequalities with cross correlations, 
which are the main result of this note. To explore to 
what extent these inequalities can be exhausted, we shall 
turn to two examples, which first called our attention to 
this problem. 

First, considering the time evolution under an inte- 
grable Jaynes Cummings Hamiltonian of a wave packet 
forming a product state, the group of M.C. Nemes Q 
observed a strong maximum in the purity after half the 
period of the classical orbit around which the packet was 
constructed 5] . It is intuitively clear that the autocorre- 
lation function will show a revival only after a full period, 
and thus cannot contribute an explanation for this phe- 
nomenon. Note that the Jaynes Cummings Hamiltonian 



H = ftwaJa + HeJ z 
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describing an oscillator with standard annihilation oper- 
ator a and SU(2) spin J, has a mirror symmetry given 
by 



—a, J± — ► — J±, J z 



J, 



(2) 



which is also evident in the figures 1 and 2 of Ref. 
We may then conjecture that the oscillations seen in the 
linear entropy, defined as one minus purity, are due to 
maxima in the autocorrelation functions for full periods 
and maxima in the cross correlation functions with the 
mirror image of the initial state for the half-period. 

Second, following the time evolution of a similar packet 
under echo dynamics with the same Hamiltonian and a 
slight detuning as perturbation, we found that coherent 
states conserved purity to higher order in h than other, 
e.g. random states Here the conjecture is that the 
cross correlation function with the classically transported 
image of the original packet will yield an explanation. 

These examples serve to illustrate our analytic results 
in the framework of a model, whose use is widespread 
in atomic physics and quantum optics, in particular for 
Rydberg atoms in cavities and for atoms in Paul traps 
with externally driven fields. The first case will become 
rather obvious once the inequalities are derived, and for 
the second case we shall derive the result within the linear 
response approximation. 

Consider a composite system with a Hilbert space 



H = H 1 ®H 2 



(3) 



consisting of two factor spaces Hi,2 which may have ei- 
ther finite or infinite dimensions. 

Let \tp) be an arbitrary pure state of a composite sys- 
tem which, after tracing over the subsystem 2, defines a 
reduced density operator over the subsystem 1 



Pi = tr 2 \yj)(ip\. 



(4) 



We can prove the following general inequality estimating 
the purity 



I[pi 



trip? 



(5) 
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of a reduced state p\ . 

Theorem: The following inequality holds 

\m)\ A <i[ P i]- (6) 

for any product state \<f>) = \fa) ® \fai)- 

Proof: Uhlmann's theorem states for pure states 
p = \if>}{ip\ and a = \4>){4>\ that 

|<<^)| 2 = trpcr < tn[tr 2 ptr 2 a] = (fa\ Pl \fa). (7) 

Then, applying Cauchy-Schwartz inequality for operators 
\tr(A^B)\ 2 < tv{AA^)tx{BB^) we get 

|tn[tr2ptr 2 <7]| 2 < tri[tr 2 /!>] 2 tn[tr2CT] 2 = I\pi]. (8) 

Combining inequalities Q and ||SJ we get the result 10. 

Corollary 1: An interesting special case of the 
above result is obtained if \fa) = | -0(0)) is an initial 
disentangled state of a unitary quantum time evolution 
\ip) = \tp(t)). Then our results says that the entangle- 
ment growth as measured by purity is bounded from 
below by the autocorrelation function |(^(0)|V'(i))| 4 - 

Corollary 2: A more general situation arises if 
\fa) is any factorizable state as required by the theorem. 
Then we obtain that the growth of entanglement as 
measured by the purity is bounded from below by all 
cross correlation functions of the type \{<f)\ip(t)) | 4 . 

The second corollary includes the first as a special 
case and is the general result we need in order to apply 
cross correlation functions to understand oscillations or 
slow growth of the entanglement. 

To explore the possibility to approach equality in the 
second corollary (and the main theorem) we have to look 
at the eigenvalues of the reduced density matrix (which 
are identical for both subspaces 0). As the purity is the 
trace of an analytic function of the reduced density ma- 
trix the optimal function <fi to establish cross correlation 
with, is a product of the eigenfunctions corresponding to 
the largest eigenvalue. Note that the eigenvalues are the 
same for the reduced density matrix in both subspaces. 
If we denote this eigenvalue by 1 — S the cross correlation 
function will take the value | (<j>\ip(t)) | 4 = (1 - S) 2 . Ac- 
tually, purity I will obey the slightly sharper inequality 

(l-6) 2 +5 2 > I > (l-5) 2 + . f 

mm(dim(/ti), dim(r£ 2 )) — 1 

(9) 

where dim indicates the dimension of the corresponding 
Hilbert space. We obtain this inequality by assuming 
the two extreme cases: either the missing intensity 6 2 
is concentrated in a single product state or it is evenly 
distributed among all such states. The former provides 
the upper bound and the latter the lower bound. As the 



lower limit is larger than the value (1 — 8) 2 we obtain 
from corollary 2 that equality in © can only be reached 
for I = 1. Note that the inequality iJBJ provides a fairly 
sharp lower bound if purity is close to 1 because then 
1 — J~ d while I — \(faip(t))\ 4 ~ S 2 , i.e. the error of the 
lower bound is quadratic in the deviation of purity from 
unity. 

We can now hope to find sets of functions fa such that 
the cross-correlation with these will explain the behavior 
of purity as long as the latter is near to 1. Yet, to be illu- 
minating we must be able to construct such a set without 
diagonalizing the density matrix. If we have a symme- 
try group which acts separately on the two subspaces, 
this set is trivially generated by applying the symmetry 
operations to the initial function "0(0) ■ For the exam- 
ple discussed in the introduction a cross-correlation with 
the function obtained by applying the reflection Q to 
-0(0), will have maxima at 1/2, 3/2, ... of the full period. 
We shall now turn to the other example we mentioned, 
namely the slow decay of purity for coherent states in 
an integrable system. This question was discussed in 
for arbitrary Gaussian packets in the context of echo dy- 
namics, and illustrated again in the Jaynes Cummings 
model. 

We can follow the line of considering cross correlations 
a little further in this context if we consider the quantity 
ti'i[tr 2 p tr 2 cr] = R{pi,a\) considered in Eq. Q. For 
p = p(t) and a — p(0) this quantity was introduced as 
reduced fidelity in [9j implying a reduced autocorrelation 
function if we do not treat an echo situation. With this 
notation we have I(pi) > R 2 {p\,a{) > |(0|0)| 4 , but 
in addition the second identity is fulfilled if cf> is chosen 
as the eigenfunction of the largest eigenvalue of p\. By 
choosing p\ = p(t) and a arbitrary, we have an additional 
cross correlation that we may use. The basic advantage 
of this quantity is, that it is, just like purity, defined on 
the subspace alone, and may therefore be handier in some 
situations. As far as approaching the identity with purity 
is concerned this quantity has no advantages. 

We shall next follow the same line of reasoning to find 
an optimized factorizable function 4> to understand the 
behavior of the purity when it is near identity for a gen- 
eral Gaussian wave-packet 

(x|V>) = Cexp(^[(x- X) ■ A(x - X) + P ■ x\), (10) 

centered at (X,P) and having the shape described by 
the (generally complex) matrix 

corresponding to a division of a d = d\ + d 2 dimen- 
sional configuration space into d\ and <i 2 dimensional 
parts x = (a; 1 ,x 2 ). The purity of a reduced wave- 
packet p{xi,Xi) = J c?x 2 (x 1 , X2\ip) (V^i: X2) is 7 = 
J dxidx' 1 \p(xi,x' 1 )\ 2 which can be evaluated in terms of 



3 



a 2d dimensional Gaussian integral 



/ = (det ImA) 



Im An 
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-1/2 



(12) 

where the vertical lines indicate the determinant of a ma- 
trix. According to our result, this expression should be 
bounded from below by cross correlation with any dis- 
entangled state which we again choose to be a Gaussian 
packet centered at the same point in phase space 



Dexp(~[($-X)-B(£-X)+P-x\), (13) 



with 



B = 



Bn 
B 22 



(14) 



Straightforward evaluation of Gaussian integral yields for 
the cross correlation 



i (det Im A) (det ImB) 



| det (A - B* 



(15) 



which, as we have shown, is strictly smaller than l|12|) for 
any B of the form l|14fl . 

One may now ask the following question: Which packet 
(|13fl optimizes the inequality 10? In general this is a 
complicated optimization problem, however we may solve 
it asymptotically when purity is close to 1. In that case, 
off diagonal elements of the shape matrix are small, so 
we may write Ai 2 = eZ\ 2l A 2 \ = eZ 2 % where e is small. 
It turns out that the inequality is optimized to leading 
order in e if we take Bn = An,B 22 = A 22 , namely in 
that case |(</#)| 4 = l-cie 2 + 0(e 4 ), / = l-c 2 e 2 + 0(e 4 ) 
where c\ — c 2 . 

We note that the cross-correlation l|15fl is manifestly 
h— independent. It may be used to compute dynamical, 



time dependent cross-correlation function of a propagat- 
ing Gaussian packet \ip(t)) (as long as linearized dynam- 
ics is a good approximation), as we know that the time 
dependence of the shape matrix A — A(t) is given by the 
ratio of two pieces of the classical monodromy matrix 
[Toj . We have thus shown that to lowest order we can 
easily optimize the reference wave function \4>) though it 
will clearly have to vary with time. We could first sim- 
ply use the transported original (and factorizable) state 
IV'(O)) centered around the new position, i.e. B = A(0), 
where {(j)\ij){t)) may be called transported autocorrelation 
function. This would not optimize the lower bound, but 
for short times it might be quite good, and it would cer- 
tainly make the bound h independent. Furthermore, one 
could apply the bound l|15|l to the case of echo dynamics 
where liner response approximation is valid even for long 
times if the perturbation is sufficiently small (see ||J). 

We may summarize our results by stating that we 
found a general inequality between the fourth power of 
cross (and auto) correlation functions and purity. Fur- 
thermore, we were able to assert, that the lower limit is 
almost identical to the value, if the purity is near to one. 
Thus any such value of purity can be explained by the 
fact that the evolving state acquires a large overlap with a 
factorizable state, and the difference between purity and 
our lower bound is of second order in the deviation of 
purity from identity. The knowledge of this factorizable 
state may contribute to the understanding of the slow 
decay or renewed maxima (also known as recoherence) 
for the purity. Particular examples of such situations are 
given in the case of symmetry group whose representa- 
tion factorizes in the variables in which we wish to split 
the spaces, and in the case of coherent states. 
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